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Abstract 

The classical trajectories for FLRW universe with varying speed 
of light are obtained for the cases in which the cosmological constant 
depends or not depend on the velocity of light. The theory is then 
quantized and the corresponding WDW equation is solved. It is shown 
that the method of causal interpretation of Bohm can be applied suc- 
cessfully to the theory. Finally the Bohmian trajectories are obtained 
and compared with the classical ones. 
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1 Introduction 

In the early works of Albrecht, Magueijo and Barrow [H [2] for presenting 
a varying speed of light theory, the lorentz invariance is broken and thus 
there is a preferred frame to formulate the physical laws. These authors 
have assumed that there is a preferred frame in which the speed of light is 
only a time dependent field and this frame is identified with the cosmological 
frame. Moreover in this theory a minimal coupling principle is assumed. 
According to it, in the preferred frame the curvature and the Einstein's tensor 
are computed by fixing the speed of light. Thus the Einstein's equations 
don't have any correction terms in this frame and one should only replace a 
fixed speed of light by a time-dependent one. Although this theory has the 
ability to solve some of the problems of the standard Big-Bang cosmological 
mo dels [H [3], it is not covariant and the comoving time is chosen as a specific 
time coordinate. 

Therefore varying speed of light theory in its minimal formulation starts 
from the Einstein-Hilbert action in which the speed of light is substituted by 
a scalar field. In addition the matter lagrangian is assumed to be explicitly 
independent of the velocity of light and a term in the lagrangian is introduced 
which depends only on the velocity of light field in order to determine its 
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dynamic. 

In order to make more simplification, Barrow [2, 4J considers cases called 
pre-set c in which the speed of light varies as a power of the expansion factor 
and so there is no dynamical term in the lagrangian for the velocity of light. 

It is shown that it is possible to generalize these ideas [5J to have a varying 
speed of light theory, in addition to preserving the general covariance and 
local lorentz invariance. The price that have to be paid for this is to introduce 
a time-like coordinate which is not necessarily equal to ct. In terms of 
and X, one has local lorentz invariance and general covariance. The physical 
time t, can only be defined when dx^/c is integrable. 

Of course introduction of a varying speed of light in the Einstein-Hilbert 
action is not unique. This is because of the fact that multiplying it by 
any power of light velocity and then introducing a varying speed of light 
would give different results. This leads a more general theory [S]. Moreover 
the cosmological constant generally depends on the velocity of light. If one 
assumes that it has field theoretical origin, it should scales as c^. In general 
one can assume c" dependence for the cosmological constant [6J. (See [7] 
and references therein for a comprehensive review of varying speed of light 
theories.) 
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Using quantum cosmology with varying speed of light, the semi classical 
tunneling probability in a minimal coupling model with pre-set c is diss- 
cused in [8] and it is shown that it is different from the corresponding one in 
quantum cosmology with constant velocity of light. Furthermore in [5] the 
quantum cosmological aspects of the minimal varying speed of light theory 
is discussed considering cosmological constant term and instantons. 

On the other hand, applying Bohmian causal interpretation of quantum 
mechanics to quantum cosmology has several positive aspects compared to 
the Copenhagen interpretation[TU | [TT | [T^ . For example, it provides time evo- 
lution of dynamical variables via Bohmian trajectories, whether the wave- 
function depends or not on the time variable. It also leads to the possibility 
to avoid the initial singularity via the quantum force. It can lead to a graceful 
exit from inflation epoch. Investiagation of causal interpretation of quantum 
cosmological models with varying speed of light could be fruitful. It is inter- 
esting to see how the admission to the velocity of light to be a function of 
cosmological time can affect the Bohmian trajectories in early universe. 

Here we shall consider the general action of [51 [7| and find it's exact 
classical cosmological solution for the fiat universe and with or without cos- 
mological constant, but without matter fields. After this we shall quantize 
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this model and obtain the exact solutions of Wheeler-Dewitt equation for flat 
and non-flat and with or without cosmological constant. We shall show that 
it is possible to present a causal interpretation of the theory. The quantum 
trajectories are calculated and compared with the classical ones. 

2 The model 

The model for varying speed of light theory we use here is the one presented 
in [5], with the action functional: 

A = J d'x,/^ (e"'^ (7^ - 2A(^) - tV^^V^V^) + e^'^£„(0„ 9^0,)) (1) 

in which = c/cq, cq is a constant velocity and we put SttG/cq = 1. We also 
have a dynamical term for the velocity of light with the coupling constant 
7, and 0j represent matter fields. In the above action one assumes that 
the c"^ factor of Einstein-Hilbert action is broken into a part behind the 
gravitational term of the action and a part behind the matter term. In 
order to remain faithful to the standard action of general relativity, one 
should set a — /3 = 4. Another point is that, here we have assumed that 
there is a time-like coordinate x", which is not equal to ct. In fact, since c is 
a field, cdt is not necessarily integrable. Therefore definition of physical time 
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is only possible when one can integrate dx^ jc [5]. 

Here we shall deal with the minisuperspace of isotropic and homogeneous 
cosmology. Denoting derivative with respect to the coordinate with a dot, 
the action functional for a FLRW model would be: 



in which V3 is the three-space volume and A(?/)), a and k are the cosmological 
constant, the scale factor and curvature parameter respectively. Depending 
on the nature of the cosmological constant one can distinguish two cases. 
First, gravitational cosmological constant which we assume that it is a 
fundamental constant of dimension {length)~'^ and hence it is ■^/'-independent. 
Second, field theoretical cosmological constant (i.e. one obtained from the 
field theoretical vacuum contribution to the cosmological constant) Aj = 
Aoe"'^. In this paper we shall choose n = 4. 
The corresponding canonical momenta are: 



A = V3 dx° (e"'^ [-6aa^ - Gaa^aip + 12ka - 2A{ip)a^ - 70^^^ 




(2) 




6a(2d + aaip)e"^ 



(3) 




2a2(7aV' + 3ad)e"^ 



(4) 
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n^, = ^ ^ 2eP^a^^ (5) 
d<i)i d(f)i 



The Hamiltonian can be calculated as: 



- Ukae""^ + 2A(V^)a2e"'^ + 2e^'f'a^T (6) 

where: 

' 12(27-3a2)2 ^'^ 

m^O^Ef'^^^-^-^l (9) 

where T is the zero-zero component of the energy-momentum tensor of the 
matter field. 



3 Classical trajectories 

Let us choose A; = in this section. We shall distinguish three cases. By a 
c-dominated universe we mean = A = 0. A (c — A^,) -dominated universe 
means that Cm = but A is not zero and is of gravitational kind, that is 
■0-independent. Finally a (c — A/) -dominated universe is one with Cm — 
and A = A/(^) = Aqc^^. 
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3.1 c— dominated universe 



The classical equations of motion in this case can be derived from the action 



via putting Cm = and A = 0. They are: 



a a 



2- + — + aV' + 2a-^ + (a^ - 7/2)^^ = 



a a' a 

3a — h 3a— + ■yip + 37--?/^ + —^ip^ = 



a 2 



These equations have the following solutions: 



(10) 
(11) 



ao 

c 
Co 



where A and -B are one of these cases: 

A = 5 = 

B 



A _ 2(a^-7) 



2a 



(12) 
(13) 



(14) 



The first solution is the flat Minkowski space-time, with constant speed of 



light. For the other two cases time can be defined as: 



dx^ (x°) 



0\1-B 



So that: 



c co(l-5) 



— = ((1 - B)cot)— ~ t^- 
ao 



(15) 



(16) 
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c , a . B B 

_ = (_)A^t— (17) 
Co ao 

Depending on the values of a and 7 these can be either decreasing or increas- 
ing functions of time. In order to solve the horizon problem of cosmology, 
one should set [7J d > and d/d — c/c > 0. This leads to some limitations 
on A and B (and therefore on a and 7): 

^ <1, (18) 



1-B ' 1-B 



3.2 (c — Ag)— dominated universe 

The equations of motion are: 



^2 



2- + ^ + ai^ + 2a-i) + {a^ - 7/2)^^ - A. = (19) 



3a- + + 7^^ + 87-^' + — - aA„ = (20) 
a a 2 

with the following solutions: 

- = eP^" (21) 
ao 

- = e^ = 6-?"° (22) 

Co 
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where p and q are one of these four possibihties: 



P 
P 
P 
P 



^6a2_9T,/2 
(a^-7)VA; 



Q 
Q 
Q 



/6A„ 



OA /A. 



(23) 



y'ea* +372-17702/2 
oa/aT 



^6o4+3^2_i7^cj2/2 ^6a4+3^2_i7^Q2/2 

Although all of these solutions are acceptable, but the third one leads to 
anti-de Sitter space-time for a = and 7 — > and the fourth one is de 
Sitter space-time in this limit. 
Time can be defined as: 



t 



_p-qx" 



qco 



(24) 



So that: 



ao 



{-qcot) 1 ~ t' 



c / a a 1 
— — — jp ~ - 

Co ao t 



(25) 
(26) 



In this case the speed of light decreases as time increases. Again the horizon 
problem can be solved provided: 



^<0 



(27) 
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3.3 (c — A/)— dominated universe 



The equations of motion are: 



... 2 

2- + ^ + at/; + 2a-4> + (a^ - j/2U^ - K^e^^ = 
a a 

3a- + 3a4 + 7^^ + 37-^^ + — - aAoe'^ - 4Aoe^^ 
a a"' a 2 

with the solution: 



Co 

where Oq and s satisfy these algebraic equations: 

Ss^ - (a + 2)s + _ = 



Qas" - 3(a + 7/2)s + 2i^±il _ + 4)Aoao = ( 
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In this case time is defined as: 

^^0 2(a;0)3/2 



t 



c 3co 
and thus: 

0, ,3 , 2s 2s 

— = (oCoi) 3 ~ i 3 

c , a _i 

— ( j 2s rv.. i 3 

Co Oo 
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Again the speed of light is a decreasing function of time, and the solution of 
horizon problem leads to: 

s > 1 (37) 



4 Quantum solutions 

The classical model of previous sections can be quantized through the method 
of canonical quantization, i.e. by setting Ha ^'^'h.i ~^ ""^^ 

— i> —i-^ in the hamiltonian ([6]). Here we set h = 1 and ignore factor 
ordering ambiguity. Since the theory is invariant under x° reparametrization, 
the Hamiltonian should be set equal to zero, and the WDW equation is then: 

da'^ a? dip'^ a dadip 
+ 2a^A(V')e2"'''<l> + 2a^e("+^)'^r(0i, -^^)^ = (38) 

in which $(a, ip, (pi) is the wavefunctional. 

Let us now try to solve the WDW equation for some cases. 



4.1 c— dominated solution 



The WDW equation would be: 



^92$ ^292$ aAa i2/taV"'^$-0 (39) 

^ da^ a? dip'^ a dadip 
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This equation can be solved via separation of variables method leading to 
the general solution: 

• k^O 

$(a,^) = e^'^^a^ (40) 

where 

A2u;^ + iA2UJas - Ais{s - 1) = (41) 
The general solution is the linear combination of this solution. 

• A; 

*(a,^) = a"V,(c2a2e"'^) or a''Y,{c2a' e''^) (42) 

where 

ci = — ^ 43 

4 = -72k (44) 

2 2Aj{2A^-a'A2 
" ^ {4A, - a^A2r ^^^^ 

4.2 (c — Ag)— dominated solution 

Here we choose k — 0. The WDW equation would be: 
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with the solution: 



where: 



-3^1 

Cl = 



a^A2 - 6^1 

2 2a^A2 - 9Ai 
^2 _ -ci(ci - 



4.3 (c — A/)— dominated solution 



Again choosing k — Q. The WDW equation would be: 
with the solution: 



$ = 0^^^(020=^6(2+")^) or a"iF^(c2a3e(2+")'^ 



where: 



3^1 - 2a 

Cl 



C2 = 



6>li - 2q;^2 - Oi^M 
2Ao 



Mi + 2(a- l)(a + 2)^ 
2 _ -ci(ci - 

~ -9^1 + 2(a - + 2)^2 
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5 The method of causal interpretation 

It is shown by de-Broglie and Bohm[T3] that it is possible to build a causal 
interpretation on the top of quantum theory. This is achieved by adding 
some axioms to quantum theory and removing some others. In de-Broglie- 
Bohm theory, the state of the system is denoted by in which q 

represents the degrees of freedom of the system, $ is the wavefunction, and 
q{t) is the trajectory of the system. The wavefunction dynamics is governed 
by Schrodinger equation, while the trajectory is defined via the guidance 
relation, dS/dq = 11. 11 is the conjugate momentum of q, and S is h times 
the phase of the wavefunction. Therefore what is added to quantum theory 
is that the state of the system is not determined by the wavefunction only, 
but one should add the trajectory of the system to it. 

On the other hand in de-Broglie-Bohm theory one does not need to 
use the projection axiom of the quantum theory. In fact, it can be shown 
easily [H] that when some property of the system is measured, the tra- 
jectory of the system is such that after a very short time the state of the 
system is one of eigenstates of the observable and the probability of finding 
any eigenvalue is given by the Born law. Therefore de-Broglie-Bohm theory 
has not the measurement axiom. 

15 



In this way, the quantum theory of de-Broghe-Bohm, which we shall call 
causal quantum theory, is really a different theory from quantum theory. 

All of these things are motivated by decomposing the Schrodinger equa- 
tion into two parts, by writing $ = Rexp{iS/h). In this way one gets two 
equations. One of them is the continuity equation with p = as the en- 
semble density and the other one is a Hamilton-Jacobi equation, with some 
additional term called quantum potential which is responsible for quantum 
behaviours. 

For our model, the causal interpretation can be derived very easily via 
putting 

$(a, i)) = R{a, ip) exp{iS{a, ip)/h) (56) 
in equation (1381) . Here we neglected matter part for simplicity. The result is: 

(57) 

n{dS/da, dS/dip; a,-^) + Q = (58) 

where the first term of the latter equation is the Hamiltonian ([6]) in which 
Ua and are replaced by dS/da and dS/dip respectively. This is the 
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Hamilton-Jacobi equation with the additional quantum potential term: 

aR \ ^ da? a^ di/j^ a dadip J 

The Bohmian trajectories can be evaluated via the guidance relations: 

f)r • (99 

n„ = — = -6a(2d + Q;aV')e°^ = — (60) 
oa da 

= ^ = -2a\^a^|; + 3ad)e"'^ = ^ (61) 
dip ow 



6 Bohmian trajectories 

Bohmian trajectories highly depends on the wavefunction of the system. For 
example if the wavefunction is real, the guidance relation leads to zero con- 
jugate momenta. But if one chooses some general linear combination of real 
solutions, one would get a non-trivial trajectory. This shows that one can- 
not speak of Bohmian trajectories and their difference with the classical ones 
without specifying that which linear combination of independent solutions of 
the WDW equation is chosen to be the wavefunction. One criteria for choos- 
ing the wavefunction could be having classical limit. That is one may fix 
the linear combination via requiring that for large times, say, the Bohmian 
trajectories tends to the classical one, in other words the quantum potential 
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is ignorable in comparison to the classical potential for large times[T^[T5]. It 
must be noted that in general, it is not possible to find solutions with classi- 
cal limit, i.e. making the quantum potential going to zero. As an example, 
it is shown recentlyfTB] that for nonseparable dynamical quantum systems 
Bohmian trajectories never can approach the classical ones. 

In what follows we shall choose specific wavefunctions and investigate 
the corresponding Bohmian trajectories for fiat universe to demonstrate how 
these trajectories can be. 

6.1 c— dominated universe 

Let us choose the wavefunction of this case as: 



(62) 



which leads to the conjugate momenta: 



n, 



■a 







(63) 



= a; 



(64) 



and thus the Bohmian trajectories are given by equations: 



a = 




a? 



(65) 
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These two equations can be combined to give the relations 



ao 

where oq is a constant of integration and 

a = -— ;t7- „ „, = C — constant (68) 

Therefore we have 

a = oo + Cx° (69) 



Co Vao + Ca;°. 
Since c is a function of x° only, the time can be defined as before: 



In terms of this physical time we have: 



(70) 



(ao + Cx°)^+'/°= (71) 



— ^ tV(l+2/a) (72) 



Oo 

^ ^ ^2/(2+") (73) 

Co 

Comparison with the classical solution shows that it could or could not be 
the classical solution. In fact there are specific choices of a and 7 that 
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the Bohmian trajectories and the classical ones coincides. This solution is a 
contracting universe which has no horizon problem (d < and a/ a — c/ c < 0). 
The quantum potential for this solution is: 

Q = -Ais{s - l)a-'e-°'^ ~ ^-VCi+s/a) (74) 
6.2 (c — Ac,)— dominated universe 

Now let us choose a more realistic wavef unction for this wavefunction 
which leads to trajectories tending to the classical ones for large scale factors. 
Using the asymptotic forms of Bessel functions: 

JJx) ~ \ — cos(x - i/7r/2 - 7r/4) (75) 

V nx 

Y^{x) ~ \ — sm{x - i/7r/2 - 7r/4) (76) 

V nx 

One observes that a linear combination of the form: 

$ = a^^ (J^(c2a3e"^) + iY^{c2ah''^)) (77) 
has correct limit. For large a^e°'^ the phase of the wavefunction is: 

S = caa^e"'^ + constant (78) 

so that: 

n„ = Scaa^e"^ (79) 
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n^, = acsa^e"^ (80) 
and the equations of motion are: 

a = ^C2{a^-j)a (81) 

ip = \ac2 (82) 



with solutions: 



Time is equal to: 



ac2 

So that in terms of time we have: 



^ = e^^2(«'-T)-° (83) 
- = e^""^^° (84) 

Co 



Cot = — ^-e-i^'^^xo ^gg^ 



Co acQC2t t 



ao V"2^°'^'J " ^'""'^ ^^^^ 
Comparison with the classical solution shows that we have the correct limit. 
The quantum potential is now: 

Q = {-A,{ci - 3/2)(ci - 5/2) - a^A2{ci - l)/2)a-3e-"'^ (88) 



(86) 
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which is ignorable in comparison to the classical potential when a^e"'^ is 



large. 

For small a^e"'^, using the expansions: 

U^) = f^^f + • • • (89) 
2^ 1 

ix) = x-"" + — cot T^vx"" H (90) 

^ ' r(l-i/)sin7ri/ r(l + i/) ^ ' 

The phase of the wavefunction is 

tan S = + a-6^e-2a.v> = ^'a-^'^e-^"'^'^ (91) 

r(l — I/) sinvri/ 

The canonical momenta are then: 

n„ = -^o^'^^-'e^'^''^ (92) 

nv^^-^a^V"^^ (93) 
and the equations of motion are: 

a = Sa^'^-'e^^^'^-i)'^ (94) 

= S'a^'^-^^e^^'"-^)^ (95) 



in which: 

C 27 - 3q;2 
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and 



C 27 - 3a2 
Dividing these two equations we have 



B'^-Z,^ (97) 



a _B ^ f c\b 



— = eW^=[-\ (98) 
ao Vco/ 

And the solution is: 

a ~ (99) 

i) ~ t-^^ (100) 

where 

7 — 0;^ 

" (37 + a - 4a2) + 2z/(4a2 - 87) 
— 

" (37 + « - 4«2) + 2z/(4«2 _ 3^) (102) 
The quantum potential can be calculated as: 

Q = (-Ai(ci-3z/)(ci-3z/-l)-aVM2-«^z^^2(ci-3z/))a-^-V"'^ (103) 

In order to illustrate how could the Bohmian trajectories be, we present here 
some plots of trajectories and quantum potential. First we choose 7 = 5 and 
a = —2. The quantum potential for this case is plotted in figure ([1]). The 
classical and quantum trajectories of the scale factor and the velocity of light 
are plotted in figures ([2]) and ([3]). It is clear that the trajectories tends to 
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classical ones for large times. In this classical limit the quantum potential 
goes to zero. 

As a second example we choose 7 = 2.5 and a = 1. The quantum 
potential, the classical and quantum trajectories of the scale factor and the 
velocity of light are plotted in figures (jlj), ([5]) and ([6]). Again the trajectories 
tends to classical ones for large times and the quantum potential goes to 
zero. 

6.3 (c — A/)— dominated universe 

Since the wavefunctions in this case are very similar to the previous case, so 
are the trajectories and we do not present it here. 

7 Conclusion and remarks 

In this paper we have obtained the Bohmian trajectories corresponding to 
a varying speed of light theory. This trajectories are shown that can have 
correct classical limits. We have done this for a matter free universe so that 
only we are able to see the status of the horizon problem for the trajectories. 
In order to disscuss about other cosmological problems like the flatness prob- 
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lem, one needs to extend the model to the universes with matter. For this we 
have to solve the WDW equation with matter fields and obtain the Bohmian 
trajectories using the method of [17]. We shall do this in a forthcomming 
paper. 

It must be noted that the specific choice of the linear combination in 
equation fl77l) imposes the boundary condition that the wavefunction behaves 
as an outgoing wave. This is what expected, because in a region where one 
has classical limit, one expects that Bohmian trajectories and the WKB path 
coincide. Therefore the appearance of an outgoing wave as the boundary 
condition is the sign of having the classical limit. 

Another important thing should be noted here that although the investi- 
gation of whether these solutions of WDW equation have unitary evolution 
needs addition of matter fields (as it is discussed in the literature, see e.g. 
|18j). but since the the linear combination of the equation fl77p have only out- 
going WKB wave limit, it would have a unitary evolution, see [IBJ. For other 
linear combinations of the solutions it is suitable to choose a wavepacket with 
compact support to have unitary evolution. 

We studied here the solutions which have the classical limit, but as the 
reader may noticed the classical singularity is not avoided e.g. in the figure 
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([2D. It is possible to choose the parameters in such a way that the classical 
singularity is avoided but the classical limit is not present. The fact that 
quantum force has the ability to stop gravity from making things singular 
is wellknown and disscused in the literature |19i UU{ [21] ■ The possibility of 
having nonsingular solutions with correct classical limit would be disscused 
in future works for more realistic cases when matter is present. 

At this end, it should be noted that we have ignored factor ordering am- 
biguity. In the causal approach to quantum cosmology, although the form 
of quantum potential depends on the regularization and factor ordering, but 
there are some general results that are independent of the specific factor 
ordering chosen. For example it is shown [22l [T71 [23] that the quantum con- 
straint algebra is independent of factor ordering. Also the general results 
obtained here do not change dramatically with respect to the factor order- 
ing. This is the reason why we have ignored it in this paper. 
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Figure 1: Plot of quantum potential as a function of the normalized scale 
factor (a/ao) and normalized velocity of light (c/cq), for a (c— Ag)-dominated 
universe. For this graph we set 7 = 5 and a = —2. 
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gure 2: Classical and quantum trajectories of the normalized scale factor 
/ ao) , for a (c — Ap)-dominated universe. For this graph we set 7 = 5 and 
= —2. Time axis is normalized to QAJ KgC^. 
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Figure 3: Classical and quantum trajectories of the normalized velocity of 
light (c/co), for a (c — A^) -dominated universe. For this graph we set 7 = 5 
and a — —2. Time axis is normalized to QAJAgCQ. 
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Figure 4: Plot of quantum potential as a function of the normalized scale 
factor {a/ Go) and normalized velocity of light {c/cq), for a (c— Ag)-dominated 
universe. For this graph we set 7 = 2.5 and a = 1. 
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Figure 5: Classical and quantum trajectories of the normalized scale factor 
(a/ao), for a (c — Ag)-dominated universe. For this graph we set 7 = 2.5 and 
a — 1. Time axis is normalized to O.SJAgCQ. 
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Figure 6: Classical and quantum trajectories of the normalized velocity of 
light (c/co), for a (c — Ag) -dominated universe. For this graph we set 7 = 2.5 
and a — 1. Time axis is normalized to O.bJAgCQ. 
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